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Abstract 

In bilinear R-parity violation (BRpV), in which the superpotential includes a 
bilinear term between the lepton doublet and the up-type Higgs superfields, a sneu- 
trino LSP can decay into pairs of heavy standard model states: W's, Z's, tops or 
Higgs bosons. These finals states can dominate over the traditionally considered 
bottom pair final state. This would lead to unique and novel supersymmetric sig- 
nals with each supersymmetric event possibly producing two pairs of these heavy 
standard model fields. We investigate this possibility and find that the branching 
ratio into heavier states dominates when the bilinear term is much smaller than the 
sneutrino vacuum expectation value for a given sneutrino flavor. When BRpV is 
the only source of neutrino masses these decays can only dominate for one of the 
sneutrino generations. Relaxing this constraint opens these channels for all three 
generations. 
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1 Introduction 



As the large hadron collider (LHC) continues to successfully probe the nature of elec- 
troweak symmetry breaking (a recent breakthrough being the discovery of a Higgs candi- 
date [1, 2]), a solution to the gauge hierarchy problem, if it exists, has evaded our efforts 
so far. An elegant candidate for such a solution and one which addresses several other 
open issues as well, e.g. dark matter and gauge coupling unification, is supersymmetry 
(SUSY). Because of its theoretical appeal, it is important to understand all the guises 
that SUSY may adopt in order to recognize it if it is produced at the LHC. 

An important open issue in SUSY, which has strong ramifications for its LHC phe- 
nomenology, as well as its cosmology, is the gauge invariance of lepton and baryon number 
violating interactions. Aside from introducing many new unknown parameters, these in- 
teractions also lead to rapid proton decay. The most common solution to address this 
issue is the imposing of a discrete symmetry, R-parity, defined as Rp = (— l)3(^-'f')+2S'^ 
[3] for a review. This forbids all the tree-level lepton and baryon number violating terms 
and also causes the lightest supersymmetric particle (LSP), which must be neutral, to be 
stable. This LSP can then play the role of dark matter and its "smoking-gun" signature 
at colliders is missing energy. However, proton decay requires the elimination of only 
the lepton number or baryon number violating terms and from a theoretical perspective, 
R-parity can be ad hoc. Furthermore, an open-mindness to possible signals at the LHC 
should push us to consider alternatives. Finally, when R-parity is violated the stringent 
constraints on superpartner masses, derived from negative collider searches for missing 
energy events at Tevatron and LHC [4, 5], can be relaxed^. 

A systematic study of all possible R-parity violating terms and their effects on phe- 
nomenology is an arduous task and furthermore, one would like a mechanism for under- 
standing why proton decay is significantly suppressed. A well-motivated solution to both 
of these issues are models which can predict the fate of R-parity. A natural framework for 
this endeavor is in the context of U{1)b-l symmetries, see [9, 10, 11] for early examples^. 
While some B — L models predict R-parity conservation [14, 15], the most minimal ones 
(in terms of particle content) require R-parity violation (RpV) [9, 11, 16, 17, 18]. Even 
some non-minimal models prefer RpV from considerations of the renormalization group 
evolution of the soft masses [10, 19]. These and other models of spontaneous R-parity 
violation, such as [20, 21], have the common feature that they can be described, in an 
effective field theory way, by bilinear R-parity violation (BRpV): the only R-parity vio- 
lating terms are the mixings between the lepton doublets and up-type Higgs doublet in 
the superpotential. This makes BRpV a powerful tool for studying possible signatures of 
spontaneous RpV. Furthermore, proton decay is highly suppressed and it is important to 
mention that a gravitino LSP can be a dark matter candidate in such models [21, 22]. 

Once R-parity is broken the LSP is no longer stable and therefore astrophysical con- 
straints on its nature do not apply [23]. Accordingly, from a purely phenomenological 

^See ref. [6, 7, 8] for details. 

^Horizontal symmetries U{l)x can be also used to construct models where the RpV couplings, arising 
from effective operators, are intrinsically small [12, 13]. 
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point of view any superpartner can be the LSP, and studies of the different possibihties 
in BRpV models (and of sleptons and sneutrinos in general models [24, 25]) and their 
relation to neutrino masses and mixings have been carried out [26, 27]. 

In this paper we extend upon previous results by considering sneutrino LSP'^ decays 
in BRpV into heavy standard model (SM) final states: W+W' , Z^Z^, h°h° and tt. 
Such states, to our knowledge, had not been considered before despite the fact that they 
can dominate sneutrino decays and can yield unique and unanticipated SUSY signals^. 
Specifically, neutrino masses force RpV to be small so its only effect is on the decay of 
the LSP. Therefore, for a sneutrino LSP, every SUSY event will eventually decay into two 
sneutrinos which could then decay into one of these heavy states. Evidence for SUSY 
might then consist of events with two pairs of Ws, Zs, Higgs bosons or tops. 

The main goal of this paper is to study the sneutrino decays into these heavy SM final 
states, which are usually due to the mixing of the sneutrino with the Higgs fields, and 
to show that they can dominate over the traditionally considered bb final state. We find 
that the latter dominate roughly when the BRpV term, is smaller than the vacuum 
expectation value (vev) of the sneutrino, Vi, for a given flavor of sneutrino, i. In the case 
when BRpV is the only source of neutrino masses, this possibility can only hold for one 
sneutrino flavor, however, when this assumption is relaxed, it can hold true for all three 
generations. Therefore, if two or more generations of sneutrinos decay via RpV, it might 
be possible to rule out BRpV as the sole generator of neutrino masses. 

The rest of the paper is organized as follows. In section 2 we discuss the generalities of 
the BRpV model, in particular those related with the neutral scalar sector. In section 3 we 
derive formulas for BRpV induced mixings. In section 4 we write the relevant couplings 
for BRpV sneutrino decays, give analytical formulas for the different partial decay widths, 
analyze the constraints on parameter space enforced by neutrino data and present our 
results. In section 5 we summarize and present our conclusions. 

2 Bilinear R-parity violation 

In what follows we will briefly describe the main features of the bilinear R-parity breaking 
model, in particular those related with the neutral scalar sector. We shall closely follow 
the notation used in [29] and assume all the parameters to be real^. Throughout the text 
matrices will be denoted in bold-face. 

In addition to the MSSM R-parity conserving superpotential (where we have sup- 
pressed relevant indices): 

W^MSSM = h^QHuU' + h^'QHj' + h^LHae' + ^iH^H^ (1) 

■^A sneutrino NLSP with a gravitino LSP would not change our phenomenological results. 
^ Gauge boson and top quark pair production through a sneutrino resonance in general RpV was 
studied in [28]. 

^This simplification does not affect our main conclusions. 
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the bilinear R-parity breaking model also contains the following terms 

WBRpV = e.^e,L^i^f , (2) 

is the SU{2) completely antisymmetric tensor, i = 1, 2, 3 runs over the SM fermion 
generations and is the R-parity and lepton number breaking bilinear parameter with 
units of mass. Consistency then requires a new set of soft SUSY breaking terms in the 
scalar potential, namely 

^BRpV = Bi ei eafs Lf . (3) 
Neglecting soft flavor mixing, the scalar potential relevant for neutral scalars is 

+ (-/i e.l\Ei - Bfie^^E^E^ + q i/f + H.c.) , (4) 

with Qz = + g'^ and = (ei, e2, €3). Electroweak symmetry is broken once the Higgs 
and slepton acquire a vev, {Ed,u) = Vd,u/y/^ and (Lj) = Vi/\/2, with v = (v^^ + + 
Si=i 2 3 ^i^)^^^ — GeV and M| = g'^v'^/A. The doublets are parameterized as 

E,= (^'] , E^=(^J] , h=(fA , (5) 



E-J ' \El 
with the neutral components given by 

E', = ^{aU Vd) , ff° = -^(a° + z^°+t;„), L, = ^{u^ + ti>^ + v,) . (6) 

Here we introduce the notation i)^'^ to differentiate the CP-even sneutrinos from the 
CP-odd. 

In the basis {S^)'^ = (o"°, (t°, z//*^)"^ the linear part of the neutral scalar potential can 
be written as 

linear ( •S'^ ) = ^ S'^ , (7) 
a=l,-- ,5 

where the t^s are the so-called tadpoles. At tree-level these are (see eqs. (4) and (7)) 
4°^ = {mjj^ + fi^)vd + Dvd - fi{Bvu + e-v) , 

tl^^ = -Bnvd + imH^+ n'^)vu- Dvu+ ^ BieiVi + \e\'^Vu , (8) 

1=1,2,3 

tf^ = Dvi + ei{~fiVd + BiVu + e-v)+ m{vi , 

where v'^ = {vi,V2,V3) and D = gzi'^d^'^'u + ^^i'^i)/^- The minimization of the potential, 
Minear = 0, rcquires the tadpoles to vanish. Thus, the vevs can be determined from the 



3 



system of equations in (8) by imposing = (a = 1, . . . , 5). In particular, considering 
only leading order BRpV terms, the sneutrino vevs can be written as 

~ -^(/iC/3 - 5j s/j) , (9) 
where m?. = +M|(c^ — 5^/2 (the tree level sneutrino mass) and = tan/3 = Vu/vd- 



3 R-parity violating mixings 

Without conserved R-parity, there are no quantum numbers to distinguish the leptons and 
sleptons from the gauginos and Higgsinos and Higgs bosons respectively. Therefore, in 
BRpV several mixings between supersymmetric and non supersymmetric particles exist: 
[i) neutralinos mix with neutrinos, {ii) charginos mix with charged leptons, {iii) Higgs 
bosons mix with the sneutrinos and {iv) charged Higss bosons mix with charged sleptons. 
These mixings are important for calculating LSP decays, especially mixings of type (i) 
as they allow to fix — via experimental neutrino data — the size of the BRpV parameters. 
Since mixings of type (iv) are not of interest for sneutrino decays, in what follows we will 
only discuss analytical approximations for mixings of type 



3.1 Neutralino-neutrino mixings 

7o £7-0 



In the basis {ip^)'^ 
written as 



with 



i\, —iXs, H^, H^,Ui) the neutral fermion mass matrix can be 



-(V^°)^MjvV'° + H.c. 



M 



N 



M 



X^" |4x3 



rT I 
Olsxs 



(10) 
(11) 



where is the neutralino mass matrix: 



( Ml -g'vdl2 g'vul2\ 

M2 gvdl2 -gVu/2 

-g'vd/2 gvd/2 -/i 

\ g'vu/2 -gVu/2 -/i / 



(12) 



and Ml and M2 are the soft masses for the bino and wino respectively. 4x3 
neutralino-neutrino mixing matrix given by 



^xi^ = ild'vi \gvi ei) 
In the Weyl mass eigenstate basis, defined as*^ 



(13) 



(14) 



^The Majorana mass eigenstates are defined as x° = [F^ F^). 
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the mass matrix becomes 

Mn = N* Mn ■ (15) 

Due to the smallness of the BRpV parameters, at order e^, Mn can be block diagonahzed 
by decomposing the diagonahzing matrix N as follows [29]: 



ui) v-r I ; \-ui i u, 

The matrix H block diagonalizes Mn to the form diag(7Vf^o, m^''^), where m'^ is the 
tree-level light neutrino mass matrix. The mixing parameters can thus be determined 
to be 

^ 2|M"^o| " 2|M^o| 

/i 4|M^o| 4|M^o| ^ ^ 

where \M^o\ = -(MiMg/i^ - 2Mi/iM^ c^^S/j - 2M2/iM^ c^^s^ tg^) (with tvK = tan6'iy, 
6w being the weak mixing angle) and 

Ai = fiVi + Vdei . (18) 

Finally the block diagonal mixing matrices Nc and in (16) diagonlize the neutralino 
and neutrino effective mass matrix, which reads 

Since this matrix has two vanishing eigenvalues it can be diagonahzed by only two rotation 
matrices, namely 

U, = U.iO^'') C7,(^lf P^) , (20) 

with 

tan^ = ^ , tan^ ^^3^^^ = . (21) 

^■^ A| " Al + A| ^ ^ 

For sneutrino decays the relevant part of the mixing turns out to be the t/J ^ block, 
that from eqs. (17) and (20) can be written as [30] 

/ -ei/fi 
C/J C = -12/ fi I (22) 
yai|A| a2|A| — es/yU a4|A| 

where ei,2 = (f^J)(i,2)jej and 63 = £3 — a3|A|/i with €3 = {Uj)zjtj. The coefficients Oj are 
given by 

""'-m^y ""'-'m^y ^^-^Mj^' ^^-^Mj^' ^^^^ 



with = g'^Mi + fif'^Ms. Taking into account eqs. (20) and (21), explicitly q {i = 1,2) 
and es are given by 

_ ei(A2 + A2) - Ai(A2e2 + A3e3) 
ei — 



£2 



|A| v/Af+Af 
A3e2 - A2e3 
v/AfTA 
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A ■ e , , 

''^W ^^^^ 

3.2 Chargino-charged lepton mass matrices and mixings 

In the bases {ip'^)'^ = (— iA^, if^^, e^^) the chargino charged lepton mass matrix is de- 
termined by the following Lagrangian 

-£,. = yf ^^U + H.c, (25) 



2 \Mc 

where = [tp^ and in the basis in which the charged lepton mass matrix is 
diagonal Mc can be written as 

V^ixl3x2 M,|3x3 ; ^ ^ 

The block diagonal matrices correspond to the MSSM chargino and charged lepton mass 
matrices whereas the off-diagonal mass matrices read 

Mrx = [-^hfv) ' = (A^^^ • (27) 

Defining the Weyl mass eigenstates as^ 

p- = Ui>- and F+ = Vi)+ , (28) 

the diagonal mass matrix Mc is obtained through the biunitary transformation 

U McV^ . (29) 

Approximate analytical expressions for the mixing matrices U , V have been discussed in 
[30, 31, 32]. We here — for completeness — describe the method. The off-diagonal block 
matrix Mrx, being proportional to the charged lepton Yukawa couplings, can be ne- 
glected, and due to the smallness of the BRpV parameters the mixing matrices U, V can 



'^The corresponding Dirac eigenstates are defined as Xi — {Fi Fp 
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Ul 0\ f I 



be written according to 

v^v...(-«;)(4. f), ,30, 

where, in first approximation in the BRpV parameters, the matrices 'Sl u block-diagonahze 
the mass matrix, namely 



(31) 



and the matrices Ul, Vr diagonahze in turn the chargino mass matrix, with the rotation 
angles given by 



From equation (31) the matrices $,l,r are found to be 

^S = M^i^iM, and Cl = (33) 

which implies is suppressed with respect $,l by a factor m^/mgusy and thus can be 
neglected (Sr = Isxs). Explicitly $,l can be written in terms of the BRpV parameters 
and the coefficients entering in the chargino mass matrix: 

= Fjg^ with \M^\ = Af^ - A.f^S2J , 

&„ = -^A.-^^. (34) 

3.3 CP-even neutral scalars mass matrices and mixings 

In the basis {S^Y" — (^d? '^ui ^i^)^ mass matrix of the CP-even neutral scalars is 
determined by the following quadratic terms 



In what follows we will discuss approximate analytical formulas for the cr° ^ 



mixing. The entries of the mass matrix in (35) involve the parameters /i, tan /3, the soft 
SUSY breaking coefficients m^. and B, and the R-parity breaking parameters e^, Vi and 
Bi. We use the minimization conditions of the scalar potential {tu4,i = 0) to remove the 
parameters B^. In doing so the matrix Mjjfj can be written in terms of the CP-odd 
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neutral scalar mass ttiao , the sneutrino masses mj>- , Mz and the BRpV parameters and 
Vi, namely 

{MhH)i2 = {MIjh)2i = -«o + Ml)c^s^ (36) 
{MHH)22 = mj,,Cp + MzSp + ^m,----^ ^ e^v,+2—— + -— ^ v,v^. 

i=l,2,3 

Where the following relations have been used 

2 2 2 I 2 I 2 2 fo^\ 

= , = "^v, + + ■^Qz'^i ■ (37) 

S2/3 o 

Note that phenomenologically consistency requires the inclusion of the one-loop correction 
in the {M]ju)22 entry, which in the limit of no stop mixing is 



2 2 

4 



See [33] for a more complete expression. 

In the absence of BRpV, Mjj^j corresponds to the neutral CP-even Higgs mass matrix 
of the R-parity conserving MSSM, which can be diagonalized via the rotation matrix 



«»h::„ :) with t.„^|;±M,,,. (39) 



The elements of the 3x3 right-lower block sneutrino mass matrix in (35) are given by 



{Mii,)(^,+2){i+2) = m% + 



I 1 



i^ii>)(i+2)(j+2) = ^glviVj + eiej (with i < j and z = 1, 2). (40) 
Finally, for the 0"° ^ — uj^ mixing 2x3 block we have 
(M|,-)i(,+2) = -/iei + M|c^"' 



5 

V 



In the mass eigenstate basis defined as 

S"^ = R^°S\ (42) 
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the Lagrangian in (35) becomes 



Cso = ^iS'YMloS"> with Mlo = R^°Mso'R^"^ , (43) 



where M'^o is diagonaL Assuming real parameters R^" can be parameterized as 

iR'r= n (^^^o^ (44) 

i=i,-,5 

where the Rij = R{6ij) are 5x5 rotation matrices. 

If the cr^u ~ ^i' mixing is small — as expected due to the smallness of the BRpV 
parameters required by neutrino data — a perturbative diagonalization of the mass matrix 
in (35) can be done. By neglecting the BRpV parameters in Mjju and -/VfJ- the rotation 
matrix reduces to 

R^ = R25 R24 R23 Ri5 Ri4 Ri3 R12 ■ (45) 

When acting on M'^o the matrix R12 diagonalizes the 2x2 block Mjjfj according to 
MfjH — diag(m^o, (where hP and are the light and heavy CP-even Higgs bosons) 
and modifies the BRpV mixing matrix Mhu- 

71^2 , f^c {Mjj-)2(i+2) + Ca(M"|j-)i(i+2)\ /.^n, 

^ U (^l..)2(.+2) - ^.(M|,,),(,+,)J • (46) 

For the mixing angle we have 612 = a. The matrices -Ri(i-|-2) eliminate the first row entries 

in (46) and up to O(e^) leave the block-diagonal matrices Mhh and Mc,u diagonal. 
The rotation angles are given by 

Sa (M'|^-)2(i+2) + Ca{Mjj 

The matrices -R2(i+2) instead eliminate the second row elements in (46) leaving again, up 
to order e^, the block-diagonal matrices diagonal. The rotation angles in this case read 

Ca{MHi>^)2{i+2) - Sa{MHu^)l{i+2) 
C^2(i+2) ^ 2 2 • y^^) 

m^o - ml. 

With these results at hand and neglecting terms of 0{6'^) the full rotation matrix in (45) 
can be written as 

50 / Ra Re 



R" ^ J) . (49) 

where and R(t account for the sneutrino and CP-even neutral Higgs components of 
the mass eigenstates, and read 

Ri> = ^^2(*^2)) ' ^ (^CQ,^^i(j+2) + Sa^2{i+2) —-50^1(1+2) " Co^2(i+2)) • (50) 
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Figure 1: BRpV induced sneutrino decay modes. In (a) and (c) sneutrino decays are 
induced by sneutrino- Higgs mixing while in figure (6) by chargino- charged lepton or 
neutralino-neutrino mixing, depending on whether the final states involves charged lep- 
ton or neutrinos. The open circles with a cross inside indicate a BRpV mixing insertion. 



4 Sneutrino decays 

With the results of section 3 we are now in a position to discuss approximate formulas for 
sneutrino decays. From now on we will consider only CP-even sneutrino decays, z/^, and 
so will drop the superscript R. Possible tree-level two-body sneutrino final states include 
fermionic modes Vi — )■ I'jl'^i ^j^k and Qk^k] electroweak gauge bosons modes z/j — )■ W^W~ 
and Z^Z^ and Higgs bosons modes i>i ^ h%°, H^H°, and H+R-. 

The goal of this paper and section is not a full study of the parameter space and 
all decays but rather to show that these heavy states can dominate the traditionally 
considered decay to bb and to identify the relevant parameter space for this dominance. 



4.1 Relevant Lagrangians 

Taking into account our approximate results for the chargino-right-handed lepton mixing, 
^ji ~ (m^/msusy), mixing with right-handed leptons is zero {^r ~ 0, see eq. (33)), the full 
Lagrangian for sneutrino decays into charge leptons is reduced to [29] 

= ^ (il^°)(.+2)i Jk Pl lk^^ - ^ C/o+2)2 K Pl I -v. + H.c. , (51) 

where the mixing matrices are given by eqs. (34) and (49). Due to the smallness of 
the BRpV induced mixing we take S'^iJ^^) ~^ ^i- Note that while the first term in (51) 
necessarily leads to two same-flavor opposite-sign charged lepton final states the second 
term can lead to different-flavor signatures. The corresponding Feynman diagrams for 
these processes are depicted in figure l((a),(6)). Being proportional to SM charged lepton 
Yukawa couplings these decays are dominated by final states involving r's. 

For sneutrino decays into neutrinos (invisible decays z/j — )• 'Ylk j^k^-'j)) the relevant 
interactions are given by 

- £(0) = hj-k {Ck,^ + C,k^) + H.C. , (52) 
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® < Vi < 

\ \ 
\ \ 
\ \ 

(a) (6) 



Figure 2: BRpV induced sneutrino decays involving Higgs final states. The open circles 
with a cross inside indicate a BRpV mixing insertion. 



where the couphng reads 

5 

Ckji = S„ A/'(j_|_4)(„_|_2) {R^ ){i+2)n, (53) 



n=l 



where Sn = (1, —1, 1, 1, 1) and 

= g' N(^k+i)i - g -/V(fc+4)2 • (54) 

The elements of the neutral fermion mixing matrix correspond to the entries of the matrix 
given in eq. (16). These interactions induce invisible sneutrino decays as the ones shown 
in figure 1(6). 

In the up and down quark mass eigenstate basis the sneutrino-quark-quark interactions 
are dictated by 

= --^/l^(il^°)(i+2)2 Uk Uk Vi - -^/if (i?^'')(,+2)l4 4 Vi ■ (55) 

As for the charged lepton final states, these decays are controlled by SM quark Yukawa 
couplings and thus are dominated by hb and tt, the last one if kinematically allowed. 
The Lagrangian for gauge boson final states is given by 

C^^""^ = gv (c/3(i?^°)(.+2)l + S^(i?^°)(.+2)2 + ^) V, h , (56) 

with = Wfj,,Z^ and gw,z = g Mw,g Mz/cw- For Higgs final states we write the 
Lagrangian involving and H^^: 

C'" = -g^,uSplvi, (57) 
where the fully symmetric coupling g(i+2)jk is given by 

gijk = ~^g'z ^ '^n {R^ ){i+2)nRjk + {R^ )jnR{i+2)k + {R^ )knR{i+2)j , (58) 
n=l 

*We do not study Di -> A" and Pi — > iJ+ decays and so do not present the Lagrangians that 
govern these interactions. 
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with Un = {ui, U2,U3. . .) = {vd, -Vu, Vi . . .) and 

5 

Rjk = (R^ )fcl ~ )j2iR^ )k2 + {R" )kn. (59) 



n=3 



The interactions in (57) for j = /c = 1,2 lead to decays of type (a) in figure 2 while for 
j = 1, A; = 2 to those shown in figure 2(6). 

4.2 Partial decay widths 

Fermionic final states are dominated by third generation quark and charged leptons. Due 
to the structure of the Higgs-sneutrino mixing, rf, hh and ti final states are possible 
independently of the sneutrino flavor, whereas f{e,fi) final states are only sizable for tau 
sneutrinos. Neglecting the final state masses, the partial decay widths for f{e,fi), rf and 
bb decays can be written as 



miG 



T{iy^ r(e, /x)) = ''^J^^ m^W.^.^ , (60) 



4^2 



TTC 



r{ui ^ rr) = — =— m^,, 
4v27rcl 



2 '"'C't^ (3,4)2 ? 



[R^ )(i+2)l — C^52'^ir 



1 2 



(61) 



m^bb) = ^^m,XR'")U2n, (62) 



where is the Fermi constant. For invisible decay modes the partial decay width, 
summing over lepton flavors, can be written according to 

k,j k,j 

with Ckji given by (53). For ti final states the phase space factors are relevant, accordingly 
the corresponding decay width reads 



For Gauge boson final states the partial decay width is given by 



(64) 



m - VV) - - 4^ f 1 - 4^ + 12^) \Ar\' , (65) 



where V = Z, W, 6z,w = 1; 2 and the amplitude Aj reads 

AY = cp{R^°)(,+2)i + s^{R^°)(,+2)2 + I . (66) 
12 



with 



and 



For V V* final states the partial decay width is given by [34] 

dr{i>i^VV*) {1-xi){1-X2) + kv{2xi + 2x2-3 + 2kv) 

— — = Ki^^vv 7-. , (67) 

dxidx2 (1 — xi — X2Y + nvlv 

Kv = ^ , Ko^vv = o \Av\ mc^ Sdy (68) 

A' 1 A' ^ 10 2 40 4 \ 

V = 1 , = ^ [u ~ Y + Y j ■ (69) 

The integration variables xi^2 lie in the ranges xi = [1 — X2 — ny, 1 — /«y/(l — X2)] and 
a;2 = [0, 1 — Ky]- The parameter •jy = Ty/My (Fy being the total decay width of 
the gauge boson V) allows a smooth transition in the threshold region where the off-shell 
gauge boson becomes on-shell (m,>^ = 2My), the calculation of F(z/j — )■ VV*) thus requires 
numerical integration over the variables xi^2 in the transition region. Outside that region, 
i.e. for rriii^ < 2My — Fy it can be written according to [35] 

where Rt{x) is a kinematic function given by 

i?T(a;) = ^ ^ arccos , — (2 - 13x + 47x^) 



y/Ax - 1 V 2x3/2 J 2x 

3 

- -(1 -6x + 4x2)logx. (71) 
Finally for CP-even Higgs bosons final states the partial decay width can be written 



as 



-> S]Sl) = ^l±^X^/\rnlrn^,ml.) U,k = 1,2) , (72) 

where A(a, b,c) = {a — b — c)^ — 46c, 5'°2 = H^, and the dimensionful coupling g(^i+2)jk 
is given in (58). For the particular case j = k {h^h^ and H^H^ final states) the width 
reduces to 



1 - 4^ . (73) 
ml 



4.3 Constraints on BRpV parameters from neutrino data 

Neutrino data plays an important role in BRpV, since BRpV parameters contribute to 
neutrino masses. Therefore, even if BRpV is not solely responsible for neutrino masses. 
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Figure 3: Numerical ranges for (i = 1,2,3) parameters and sneutrino vevs as required 
for explaining neutrino data for the R-parity conserving parameters discussed in the text. 
The horizontal (vertical) solid lines in the right-panel plot indicate the values where A3 is 
dominated by V3 (t^), i.e. where A3 ~ /if 3 (A3 ^ vco^Pt-^). 



e.g. [11, 21, 37, 38, 39], one must not saturate neutrino masses via BRpV. Therefore it 
would be worthwhile to discuss this correlation briefly here. We proceed by assuming 
BRpV as the sole source of neutrino masses and later consider relaxing this assumption. 

At tree-level, BRpV allows for only one massive neutrino, as mentioned earlier. One- 
loop contributions then leave only one generation massless. Furthermore, for the approx- 
imations made in this paper, the tree-level mass is larger than the one-loop mass. The 
upshot of all this is that only allow so-called normal hierarchy neutrino mass spectrum is 
allowed here. Recent neutrino data [40, 41, 42] then allows the following neutrino masses: 



rrn =0, ma = y Am^^ = 0.00873 eV, = ^ Am^^ = 0.0505 eV. (74) 

In terms of the model parameters 

The experimental values of the atmospheric and reactor angles yield two additional con- 
straints given by eq. (21). Therefore, the atmospheric sector entirely fixes the Aj param- 
eters. 

Constraints on the e parameters arise from the solar sector. As long as the one-loop 
contribution to neutrino masses is smaller than the tree-level one (an assumption we use 
throughout), Am^^^"^ is accurately determined hj b — b one-loop diagrams [30]: 
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where 6^ stands for sbottom mixing, ei,2 are defined in eq. (24) and ABq^^^ = Bo{0, m?^,ml)- 
i?o(0,m? ,ml) (with Bo{0,x,y) a scalar Passarino-Veltman function [43]). Note that due 
to nib < "^6,2, 



Moreover the solar mixing angle can be written as [30] 



e2 



tan2C"'' = ^. {n 



Thus, eqs. (76) and (78) provide two constraints and determine, from eqs. (24), ei,2 (ei,3 
or 62,3) as a function of €3 (62 or ei). Once the A's and e's are fixed the sneutrino vevs are 
automatically fixed as well through 

Ai - Vdei 

Vi = , (79) 

see eq. (18). This in turn fixes all the relevant mixings for sneutrino decays (x°— z/, x~"~^L) 
0"^ „ — u^) once the R-parity conserving supersymmetric parameters are specified. With 
this knowledge in hand, we start to explore the consequences of the neutrino sector on 
the BRpV parameters. We note that our results have been verified using SPheno [44]^. 

Figure 3 shows typical values for ei,2 and sneutrino vevs as a function of €3. The plots 
were obtained by fixing 9^ = vr/lG, tan/? = 10 and 

(/X, Ml, M2, m^^, m^^, TUA, m^^ = (650, 550, 600, 10^ 700, 10^ 300) GeV, 

where we have assumed a conservative lower bound on the sneutrino mass of 100 Gev, 
see [45] for more details. We will use this parameter point throughout the paper. Neutrino 
observable (Am32, Am2i, %) are varied in their 3a experimental range [40, 41, 42]. We 
further assume that €3 is positive. For these assumptions we see from the left-hand 
side that ei and €2 have a lower bound and from the the right-hande side that this lower 
bound is larger than vi and f 2 respectively. On the other hand, since €3 is undetermined, f 3 
dominates it up to about |e3/yu| = 10~^. While allowing negative values for the parameters 
changes this picture, it leaves one important qualitative property the same: only in one 
generation, j, can the sneutrino vev, Vj, be larger than the bilinear mixing 
term, Cj, when BRpV is the sole contributor to neutrino masses. This will have 
important consequences in the next section. 

Also of note in the left-panel plot, €3 obeys only an upper bound determined by 
the condition {m'^y'^'^'^ > (m^'^)^"'°°P, which in terms of the bilinear R-parity violating 
parameters translates into | A| > |ep. In contrast the ei_2 parameters, due to solar neutrino 
physics constraints, are forced to lie in a "narrow" range and are such that a region where 



^We thank Werner Porod for fixing a bug in the decay routines which allowed us to calculate sneutrino 
decays to gauge bosons, discussed in section 4.4 



15 




Figure 4: Values of R-parity breaking parameters for which BRp V is not the sole source 
of neutrino physics. The parameters are such that Am^^^^ fits the experimental 3a range 
but Am2^^^ falls below the measured value (see the text for more details). 



ei,2 ^ £3 exist. Consequently, while Ai_2 are mostly determined by ei,2, A3 is controlled 
by V3 in the region where ea/yU ^ 10~^, as demonstrated by the horizontal solid lines in 
figure 3 (right-panel) which correspond to A3 ~ fivs. 

The bilinear R-parity breaking parameters selected as described above satisfy neu- 
trino data, and thus lead to BRpV models that can account for neutrino masses and 
mixings. However, it might be that these parameters are not sufficiently large to account 
for the neutrino mass scales. In that case their contribution to the atmospheric and so- 
lar masses are still determined by eqs. (75) and (76) but are such that, for example, 
Arri^^^^ < {<^)Am^^'^ and Am^^^"^ < (^)Am2f . Figure 4 shows the results for an 
illustrative case where the BRpV model fits the atmospheric mass scale as well as the 
atmospheric and reactor angles in their 3a experimental range, but the contribution to 
Am2i is subdominant^*^. Note that we plot only ^1^3, as we have found that V2 ~ V3. This 
result as well as vi <^ in the region of small 61^3 are due to the constraints arising from 
fitting 623 and ^^13. 

In summary, if data is not explained by the BRpV parameters all the e's can be small 
mainly due to the absence of the solar data constraint, thus implying that, in this case, a 
region where the three sneutrino vevs are large (Aj ~ /it>j) exist. 

4.4 LSP sneutrino phenomenology 

We are finally ready to address the main aim of this paper: showing that sneutrino LSP 
decays into heavy SM final states, i.e. W~^W~ , ZZ, h^h^, and ti, can dominate over the 
traditionally considered lighter states. Of the latter, bb has been shown to generally have 

^°Models of this kind have been considered in the hterature, see e.g. [36, 37]. 
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Figure 5: Results of a scan over SUSY parameters, as shown in the text, comparing the 
sneutrino decay width into W^W^ and bb final states left plotted in the ei-Vi plane. It has 
been assumed that BRpV does not saturate neutrino masses. Blue (red) points indicate 
where the W^W~ (bb) final state dominates. This plot therefore shows an important 
general result: the W~^W~ dominates only when the sneutrino vev is much larger than its 
bilinear mixing parameter. 



the largest partial width [25] and we compare all partial widths to it. 

We begin with a very general study ignoring all neutrino constrains (SUSY parameters 
can always be chosen in such a way so that BRpV does not saturate the neutrino masses) 
and exploring in the ej-f, parameter space where W^W~ final states dominate the bb 
decays of i>i. Without neutrino constraints, the decay properties are independent of 
generation. We scan over the following values: 

tan /3 = 2 - 50, = 500 - 1000 GeV, m^^ = 161 - 500 GeV, 
Vi = IQ-^ - IQ-* GeV, ei = 10'^ - 10'^ GeV, 

and /i below 1000 GeV but large enough so that the sneutrino is the LSP. We also use 
rrih = 125 GeV as suggested by recent LHC results. The the results are displayed in 
figure 5 where blue dots indicate the points at which the W^W~ width is much larger 
than the bb width, while the red points show the opposite. It is striking that the ratio 
of the partial widths are relatively independent of R-parity conserving SUSY parameters 
and depend only on the BRpV parameters. From this figure, one can conclude that the 
W^W~ partial width dominates the bb width only when the sneutrino vev is much larger 
than the bilinear mixing parameter and, as will be shown below, this condition holds for 
the other heavy SM final states. Couple this with the main result of the previous section: 
satisfying neutrino masses solely through BRpV means that Vi ^ can only be satisfied 
in one generation, indicates that heavy final states can only dominate the sneutrino BRpV 
decays for one generation of sneutrinos. Of course, this generation must also be the LSP 
in order for these decays to be observable. 
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Figure 6: Partial decay widths for z>i (left) and (right) to bb, h^h^, WW and ZZ 
final states as a function of lea/yu]. BRpV parameters have been fixed by using neutrino 
oscillation data (see the text for more details). Decays o/z>2 are not shown since they are 
very similar to those of Vi. 

We continue by considering the specific SUSY point discussed in the last section and 
examining sneutrino decays into 66, W~^W~ , ZZ and h^h^ assuming that BRpV is solely 
responsible for neutrino masses (the tt channel is suppressed due to the off-shell top for 
this sneutrino mass). The results are displayed in figure 6 versus e^, where it can be seen 
that for z>i, the dominant decay mode always corresponds to bb, the h^h^ decay branching 
ratio barely reaches values of ~ 10~^ (this holds for z/2, which was not plotted since its 
features are similar to ui). For the situation, however, is quite different. While in the 
region of large the bb mode dominates, in the region of small es, h^h^ ZZ and W^W~ 
final states become the dominant channels, eventually exceeding the bb mode by almost 
two orders of magnitude. 

Figure 6 supports the argument made in reference to figure 5. Reiterating: for a 
sufficiently heavy sneutrino LSP, even if bilinear R-parity violation accounts for neutrino 
data, the bb channel is not necessarily the dominant decay mode. Therefore, observing 
a sneutrino decaying dominantly to h'^h^ (or W^VT)^^ does not rule out bilinear R-parity 
violation as responsible for neutrino masses and mixings. However, if the sneutrino mass 
splittings do not allow the heavier sneutrinos to decay into the LSP, two or all three 
sneutrino generations will decay via R-parity violating couplings. In this case observing 
at least two different sneutrino fiavors^'^ decaying dominantly into heavy SM final states 
will prove that bilinear R-parity violation does not account for neutrino data, because 
this would requires Vi ^ for those generations, which is not possible if neutrino physics 

^^Although final states have a larger decay branching fraction, depending on the experimental 

environment, WW modes could be an easier experimental target, as it turns out to be the case at LHC. 

^^In principle for sneutrinos stemming from chargino decays the corresponding sneutrino flavor could 
be tagged [25]. 
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Figure 7: Decay branching fractions for final states with "strong" kinematical thresholds 
normalized to BR{i)^ — )■ hh) as a function of mj^.^ . The Branching ratios have been calcu- 
lated in the small region ({e^/fil = 10~^) and for bilinear R-parity breaking parameters 
fixed by neutrino parameters according to their best fit point values [40, 4h 4^]- 



is determined solely by BRpV, see figure 4. 

Focusing further on the heavy SM final states of 1)3, we study its decays as a function 
of its mass in figure 7 allowing a heavy enough z/3 so that even the tt channel is opened. 
All the R-parity conserving parameters have been fixed according to the supersymmetric 
point used in the previous section, the R-parity breaking parameters have been fixed by 
adjusting neutrino observables to their best fit point values [40, 41, 42] and es = /i x 10^^. 
We have checked the ratios are quite insensitive to changes in tan/3 in the range [2,30]. 
Calculations of the WW* and ZZ* modes are performed to show where these overtake 
the bb channel, close to the relative thresholds. However, the off-shell calculations for 
the Higgs and top states have been neglected since in that region the gauge boson states 
dominate. 

Note that this analysis has been done by fixing the bilinear R-parity breaking param- 
eters via neutrino data. If one sticks to BRpV models that do not account for neutrino 
masses (as in the second case discussed in section 4.3) the results will not drastically 
change, as what determines the relative importance of the ratios BRj^j^/BR^^ is the con- 
dition Vi ^ ej. However, in that case, a quantity that becomes relevant is the sneutrino 
decay length L{i>i) since it can be that the size of the bilinear R-parity breaking parame- 
ters leads to a sneutrino decaying out of the detector. Neglecting the Lorentz boost factor 
and considering only the leading processes (the ones discussed in figure 7), we have found 
that as long as the R-parity violating parameters fit data, in general L{i)i) is well below 
1 mm (inline with [27]). We have also checked that if r^, = Am-gg^lf'/Amg^gi = 10"'^ the 
decay length is generically below ~ 10 cm. Note that values of in the range (10~^, 1) 
might be in conflict with neutrino data when the contributions from the mechanism re- 
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Figure 8: Estimated upper limit for the sneutrino decay length for different choices of the 
bilinear R-parity breaking parameters. The BRpV contribution to the atmospheric and 
solar mass scales is specified by the parameter = / Am^2$i, where Amg^ai 

have been taken according to their best fit point values [40, 4I, 4^]- 



sponsible for neutrino masses are taken into account (for example a standard seesaw). 
Once the contributions to the atmospheric and solar mass scales fall below = 10~^ 
the decay length starts exceeding 1 m for certain sneutrino mass ranges. Figure 8 shows 
the upper bound on the decay length (we consider contributions only from bb and the 
heavy SM states) for = 1, 10~^, 10~^, 10^^. The downward spike at m^g = 125 GeV is 
due to the singularity of the neutral scalar mixing at m^;. = (see eq. (48)) while the 
remaining irregularities are a consequence of the different kinematical thresholds {WW, 
ZZ, h^h° and tt). 



5 Conclusions 

BRpV is a good effective theory for models of spontaneous R-parity violation, in which 
proton decay as well as the number of new parameters is under control compared to 
explicit R-parity violation. In the presence of BRpV, a sneutrino LSP decays in some 
way resemble a heavy Higgs and include light SM states, as well as heavy states which 
have not been considered in the literature before. We have studied these decays into 
W^W~, ZZ, h^h^ and ti and found that long as ^ Vi, the heavy SM modes dominate. 
As discussed in section 4.3, for models where the BRpV parameters are fixed by neutrino 
data, due to the constraints from the solar sector, this is possible only for a single sneutrino 
flavor. Therefore, if 2 or more sneutrino generations are degenerate enough to decay via 
BRpV into these heavy SM states one could rule out BRpV as the sole source of neutrino 
masses. 

For models where the bilinear R-parity breaking parameters do not contribute signif- 
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icantly to neutrino masses, the constraint <C Vi for all flavors is viable. Accordingly, in 
these BRpV models large sneutrino LSP branching fractions to l^'^VF", ZZ, hPhP and 
ti are possible, regardless of the sneutrinos mass spectrum. In this case, however, special 
attention has to be paid to the sneutrino decay length L{i)i). We calculated upper limits 
for L(ui), finding that as long as the BRpV contributions to neutrino masses are not below 
~ 0.1% the decay length is genericaly below ~ 10 cm. 

Regardless of these considerations, sneutrino decays into heavy SM final states 
ZZ, h^h^ and tt) are an interesting phenomcnological possibility that has not received 
much attention before. Because BRpV only affects the decays of the LSP, assuming a 
sneutrino LSP means that SUSY events could cascade decay to two pairs of these heavy 
SM states, which constitutes a unique and novel SUSY signature. 
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